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ABSTRACT 

We present the study of the CoRoT transiting planet candidate 101186644, also named LRc01_El_4780. Analysis of the CoRoT 
lightcurve and the HARPS spectroscopic follow-up observations of this faint (my = 16) candidate revealed an eclipsing binary 
composed of a late F-type primary (Teff = 6090 + 200 K) and a low-mass, dense late M-dwarf secondary on an eccentric (e = 0.4) 
orbit with a period of ~ 20.7 days. The M-dwarf has a mass of 0.096 + 0.01 1 M©, and a radius of 0.104^^ ^^^ R©, which possibly makes 
it the smallest and densest late M-dwarf reported so far. Unlike the claim that theoretical models predict radii that are 5% - 15% 
smaller than measured for low-mass stars, this one seems to have a radius that is consistent and might even be below the radius 
predicted by theoretical models. 

Key words. Planetary systems - Stars: individual: CoRoT 101186644 - binaries: eclipsing - Techniques: photometric - Techniques: 
radial velocities 

1, Introduction Eclipsing binaries (=EBs) with low mass ratios, like M 

.. ,^ .. dwarfs that eclipse A-K dwarfs, draw p articular attention (e.g., 

CoRoT IS the firs t space-miss ioHto search for tra ns^^Mig |Bentlev et allUM. iTriaud et al.ll2oM . The high-quahty pho- 
extrasolar p anets dBag^^^ M l Auvergne et alJ M tometric and spectroscopic data p roduced for these systems 
| Deleuil etal. ||Mll|)-^Since its launch in December 2006, 26 tran- by photometric surveys like HAT (Bak os et al.ll200l. W ASP 



sitmg planets have been discovered and many other cand idates jpoUacco et al.ll200 6). CoRoT, Kepler (B orucki et all2010h . and 

await confirmation fetzold et al. | 201 2t ^ .wa et al. | 201 2|) . The Lir follow-uiprograms, allow for precise deteniinaZ of 

discovery process of a new transiting planet includes careful j^^j properties, such as the masses and the radii of both 

^ analysis of the lightcurve (Caipano et aL 2009), as well as pho- components. This in turn enables addressing fundamental prob- 

H • t^t nc and spectroscop i c follow-up observa hons dDeegetalJ ^ ^^^j^ j^.^^^^ ^^^^^ ^^^j^^. ^^^^ ^^^^ 

_cp_-^^^^Mm^E^^Mmmm. oredcal models predict 5% - 15% s maller radii than measured 

Lightcurves with periodic transit-hke signals caused by foHgw^ss stars (e.gJTorres et al. 2010: Morales et al...201Q: 

phenomena other than a transifing planet are usually referred [Rraus et alJl201lL and references therein), 
to as false po sitiv es or false alarms (e.g., Brown 2003t 

ISanteme et alJ 120121) . In the CoRoT sample, for instance, the This paper presents the study of t he CoRo T transiting-planet 
main source of such alarms are eclipsing b inary systems in var- candidate CoRoT 101 186644 (CabreraetaD|200£), also named 
ious configurations (lAlmenara et al.ll2009h . Recently, however, LRc01_El_4780. The CoRoT lightcurve and the HARPS spec- 
such systems are becoming objects of interest in and of them- troscopic follow-up observations have led to the conclusion that 
selves (e.g., Pont et al. 2005a; Tal-Or et al. 201 1). LRc01_El_4780 (= C4780) is an eclipsing binary composed of 

a late E-type primary and a low-mass dense M-dwarf secondary 

* Based on observations made with the 1-m telescope at the Wise star on an eccentric (e = 0.4) orbit with a period of ~ 20.7-days. 
Observatory, Israel, the Swiss 1.2-m Leonhard Euler telescope at La- 

Silla Observatory, Chile, the IAC-80 telescope at the Observatory del Section 2 gives some details about the star and presents the 

Teide, Canarias, Spain, and the 3.6-m telescope at La-Silla Observatory CoRoT lightcurve. Section 3 describes the ground-based follow- 

(ESO), Chile (program 184.C-0639). up observations we performed. Section 4 presents our analysis 
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-dwarf on an eccentric orbit around an F-star 



Tal-Or et aL: Transiting low-mass dense M 

Table 1. Coordinates and magnitudes of C4780 and of the main 
contaminator in its CoRoT-photometric mask. 



LRc01_El_4780 


CoRoT ID 


2MASS ID 




101186644 


19265907 + 0029061 




RA (J2000) 


Dec (J2000) 




19^^26'"59^08 


00°29'06".4 




Filter 


Magnitude 


Source 


B 


17.09 + 0.15 


ExoDaf 


V 


16.05 + 0.07 


ExoDat 


r' 


15.67 + 0.14 


ExoDat 


V 


14.93 + 0.03 


ExoDat 


J 


14.03 + 0.03 


2MASS^ 


H 


13.63 + 0.02 


2MASS 


Ks 


13.56 + 0.03 


2MASS 


Main contaminator 


CoRoT ID 


2MASS ID 




101186975 


19265918 + 0029008 




RA (J2000) 


Dec (J2000) 




19^26'"59^21 


00°29'01".0 




Filter 


Magnitude 


Source 


V 


19.08 + 0.37 


ExoDat 


r' 


18.57 + 0.27 


ExoDat 



^Deleuil et al. (2009) 
^ Cutrietal. (2003) 
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Fig. 1. C4780 lightcurve normalized by its median flux value. 
Top: the original lightcurve. Bottom: the detrended lightcurve 
moved down by 0.05, for clarity. 



of the available data to derive the main physical parameters of 
the system. The astrophysics of the secondary is discussed in 
Section 5 within the context of the M-dwarf radius problem. 
Finally, Section 6 presents some more general conclusions. 



2. The CoRoT lightcurve 

C4780 was observed by CoRoT during the LRcOl run from May 
19, 2007 to Oct. 12, 2007. Table[I]Hsts some basic information 
on its stellar properties. Since it is a relatively faint star, only 
monochromatic photometric data were recorded. Eight transit- 
like events with a depth of ~ 0.86% an d a period of 20 684 
days were detected in the lightcurve dCabrera et al.l l2009l) . It 
was therefore identified as an interesting candidate for a Jovian 
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Fig. 2. The 8 transit-like events detected by CoRoT in C4780, 
each moved to its own timeframe (i.e. (n • P + Tq) days were 
subtracted from the time stamp of each point, where n is the 
transit number, Tq = 2454240.3 147 BJD, and P = 20.6837). 
Subsequent events were moved down by 0.02 for clarity. 

plane t residing in the s o-called "period valley" (e.g. Jjones et al.l 
l2QQ3l:IUdrvetaDl2003h . 

Figure [T] shows the original and detrended lightcurves of 
C4780, both normalized by their median flux value. Figure |2] 
zooms on the eight transit-like events detected by CoRoT in the 
lightcurve, each shifted to its own timeframe. Stellar and sys- 
tematic variability were removed from each transit by fitting a 
third-degree polynomial to the out-of-transit points. 

3. Ground-based follow-up observations 

As for other CoRoT transiting-planet candidates, we performed 
a sequence of follow-up observations to understand the true na- 
ture of C4780. Photometric observations were done to verify that 
the transits indeed occur on the main star inside the large-area 
photometric mask of CoRoT. Using spectroscopic observations 
we tried to distinguish it from other possible false-alarm sce- 
narios. Together with the CoRoT lightcurve, these observations 
were used in determining the system's nature as an eclipsing bi- 
nary system and in deriving its main physical parameters. 

3.1. Photometric observations 

Figure [3] shows images of the sky around C4780, taken by 
CoRoT and by the 1-m Wise-Observatory telescope, Israel. It 
can be seen that, in addition to C4780, the CoRoT-photometric 
mask contains at least three fainter stars, but only two of them 
are bright enough and fully contained in the mask to be able 
to cause the ~ 1% transit-like signals — the two stars ~ 6'' 
south-southeast of C4780. The brighter of the two (CoRoT ID 
101186975) is fainter than C4780 by 2.9 ± 0.3 mag in the r'- 
filter (see Table[T]), while its close neighbor is fainter than C4780 
by 4.2 ± 0.5 mag in the R-filter, as measured using the Wise- 
Observatory data. The contamination level inside the mask is 
thus -10%. 

ON-OFF observations of C4780, in which short timeseries 
during a transit and outside of it are observed and compared pho- 
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Fig. 3. Images of the sky around C4780 (the star is at the center). 
Right: R-filter image taken on August 9, 2010, by the 1-m Wise- 
Observatory telescope, Israel. Left: Image taken by CoRoT, at 
the same scale and orientation. The rectangular line at the center 
outlines the CoRoT-photometric mask. Positions of nearby stars 
are indicated by small crosses. 

Table 2. HARPS spectra and RVs (in km s'^) of C4780. 



HJD 


RV 


error 


Texp (sec) 


S/N 


2454638.741570 


16.4228 


0.0896 


3600 


2.7 


2454643.861050 


14.0218 


0.1637 


3600 


3.1 


2455304.897570 


14.3876 


0.0921 


3600 


4.7 


2455338.784230 


18.9455 


0.1095 


3600 


3.4 


2455351.719640 


27.1110 


0.1004 


3600 


3.8 


2455352.737020 


28.0719 


0.0737 


3600 


4.4 


2455353.891020 


26.5536 


0.0510 


3600 


6.3 


2455389.722540 


15.6913 


0.0853 


2700 


3.9 


2455397.769130 


22.2351 


0.1971 


2700 


2.0 



tometrically (iDeeg et al. were performed during August 

2010 using the 1-m telescope at the Wise Observatory, Israel, 
and the Swiss 1.2-m Leonhard Euler Telescope at La-Silla 
Observatory, Chile. Both observations indicated that the transits 
occur on the main target in the mask — C4780 - and not on one 
of the two contaminators. 

In an effort to recover its ephemeris, C4780 was observed 
again on May 27, 2011, using the IAC-80 telescope at the 
Observatory del Teide, Canarias, Spain, and the Swiss 1.2-m 
Leonhard Euler telescope at La-Silla Observatory, Chile. A clear 
ingress was detected in these observations, which permitted the 
refinement of the period estimation to 20.68369 ± 0.0001 1 days. 

3.2. HARPS spectra and radial velocities 

Spectroscopic observ ations of C4780 w ere done with the 
HARPS spectrograph (iMavor et al.l l2003l) mounted at the 3.6- 
m ESO telescope, Chile, as part of the ESO large program 
184.C-0639. HARPS was used with the observing mode obj_AB, 
in which the background- sky spectra were recorded using the 
second fiber (fiber B). Overall nine spectra were recorded be- 
tween June 21, 2008 and July 20, 2010 with exposure times of 
45 -60min. T he HARPS data were reduced with the online stan- 
dard pipeline (iBaranne et al.l 1 9961: IPepe et al.ll2002 l). The signal- 
to-noise ratio (= S/N) per pixel at 550 nm is in the range of 
2-6.3, since C4780 is on the faint end in magnitude for HARPS. 
Radial velocities were obtained by performing weighted cross- 
correlation with a numerical G2 mask. The derived radial veloc- 
ities (= RVs) are given in Tabled 
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Fig. 4. Three parts of the co-added HARPS spectrum of C4780 
(blue), together with the fitted model spectrum (black). Top, 
middle, and bottom panels show 32 A long parts of the spectrum, 
focused on the H a, Mg lb, and HJS lines, respectively. 

Table 3. Parameters from ground-based observations. 



Parameter 


Value 


Photometric follow-up 


Orbital period (P) [day] 


20.68369 + 0.00011 


Time of center of primary transit (Tq) [BJD] 


2455708.867 + 0.008 


Spectrum modeling 


Effective temperature (T^eff) [K] 


6090 + 200 


Surface gravity (log g) [cgs] 


4.4 + 0.2 


Metallicity ([Fe/H]) [dex] 


+0.2 + 0.2 


Projected rotational velocity (v sin /) [km s~^] 


3 + 2 


Mass (Ml) [Mo] 


1.2 + 0.2 


Age* [Gyr] 


< 7 


Broad-band photometry modeling 


Effective temperature (reff) [K] 


5800 + 400 


Surface gravity (log g) [cgs] 


4.45_o,8^^-^^ 


Distance (d) [pc] 


1100_i5o"^'' 


Extinction (Ay) [mag] 


1.1-0.4^^-^ 



* 2 cr upper limit 



The spectra were a lso analyzed with TODMOR (e.g., 
IZucker & IMazet] 119941: IZucker et al.l l2003l l2004b . a two- 
dimensional correlation algorithm customized for detecting 
faint- secondary companion in a spectrum. However, the signa- 
ture of neither the secondary in C4780 (= C4780 B) nor any other 
stellar contaminant were detected. 

To determine the fundamental atmospheric parameters of 
the primar y in C478Q ( =C478 QA), we used the method de- 
scribed by iBruntt et aP (l2010bh . which has become the stan- 
dard method for characterizing CoRoT targets since CoRoT-3b 
(e.g.. iDeleuil et al.ll2Q08l: iFridlund et al.ll2QlQh . In this method 
the observed spectra are first co-added to create a single master 
spectrum, and then a synthesized spectrum is fitted to this mas- 
ter spectrum using either the Spectrosco py IMade Easy (SIME, 
Valenti & Piskunovl [T996l: IValenti & Fisc her 2005) or the VWA 
(Bruntt et al. 2008) spectral analysis packages. 

The atmospheric parameters found this way are listed in 
Table [3l The relatively large uncertainties are a result of the low 
S/N of the spectra and the addition of possible systematic errors 
(see e.g., Bruntt etal. 2010a, 2012; Torres et al. 2012). Figure 
|4] shows three parts of the master spectrum of C4780, 32 A long 
each, focused on the H a, HyS, and Mg I b lines, together with the 
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Fig. 5. Two sets of stellar isochrones from I Pemarque et al.l 
(I2QQ4 . one for [Fe/H] = 0.05 (solid lines) and one for [Fe/H] = 
0.39 (dashed lines). Both sets are for [a/Fe] = and ages of 
0.2, 0.4, 1, 2, 4, 8, and 10 Gyr (going from left to right along 
the logg = 4.3 line). The estimated T^s and logg of C4780A 
are marked by a star and their uncertainties are marked by an 
ellipse. 

fitted model spectrum. The extended wings of B aimer lines can 
be used to constrain Tqq, while log g of late-type sta rs can be de- 
term ined from pressure-sensitive lines like Mg I b (iBruntt et al.l 
l2Qr0bl) . 

4. System parameters 

4.1. Isochrone fitting 

The primary mass and the system's age were estimated using the 
atmospheric parameters derived from the HARPS spectra and 
a grid of stellar isochrones (lYi et al.l [200 U iDemarque et al.l 
|2004|) . This was done by taking into account all age and mass 
values that fall into the ellipsoid in the (T^eff , log g, [Fe/H]) space 
defined by the atmospheric parameters and their errors. To illus- 
trate the process Figure |5]shows two sets of Y^ stellar isochrones 
of 0.2- 10 Gyr, one for [Fe/H] = 0.05 and one for [Fe/H] = 0.39. 
The ellipse defined by the estimated Tq^, log g, and their uncer- 
tainties is also shown. A lower limit of 0.2 Gyr on the age was 
set to ignore possible pre-main sequence solutions. This proce- 
dure yielded a mass estimate of 1.21 ± O.IOM© and an upper 
limit of ~ 7 Gyr for th e system's age (at a 2 cr confidence level). 
Following lBasu et aH C2Q12) we have conservatively doubled the 
mass errors to take possible uncertainties in stellar model param- 
eters into account. 

The consistency between the atmospheric parameters de- 
rived from the HARPS spectra and the broad-band photome- 
try (listed in Table [TJ was also checked. This was done by fit- 
ting the distance and extinction to minimize the between 
each isochrone point, which also contains predictions for the 
true My and color values, and the observed magnitudes. ExoDat 
and 2MASS magnitudes were translated to the systems used in 
the isochrones using the relations given bv Deleuil et al.l (l2009l) 
and Carpenter (2001). The average extinction law (Ry = 3.1) 
of i Savage & Mathis (1979) were assumed, together with the 
ICardelH . Clayton, & Mathis (1989) total-to- selective extinction 
ratios. The parameters were derived taking all model-points with 

^ xiiin + 1 i^to account. Only isochrones with metallicity in 



the range allowed by the spectrum modeling (-hO.2 ± 0.2) were 
considered. The results, given at the bottom of Table [3] are con- 
sistent with the parameters derived from the HARPS spectra. 

4.2. MCMC analysis 

The rest of the system parameters were derived by modeling the 
CoRoT lightcurve and the HA RPS RV data simu ltaneously, us- 
ing MCMC analysis (e.g., Teg mark et aP |2004 App. A). The 
MCMC run consisted of 10^ accepted steps. The input data were 
the eight-transit data shown in Figure |2l seven one-day-long 
parts of the detrended lightcurve where the secondary eclipses 
might have occurred, and the nine HARPS RVs. Errors for the 
photometric data were derived from the scatter of the out-of- 
transit points, which was found to be ~ 0.25%. 

The model consisted of 12 free parameters — P, To, r^, 
k, L3, X, Up, e cos oj, e sin oj, K, and 7, all detailed in Table IH A 
lower bound of 0.07 was set on L3 to fit the estimated contami- 
nation value (see Section 3.1), and Up was constrai ned to bebe- 
tween 0.56 and 0.66 according to the values given in [Sh^(l2QTQh 
and the atmospheric parameters derived from the HARPS spec- 
tra. Besides these two limitations, flat prior distributions were 
used. 

At each MCMC step, the light curve model was calculated 
using EBOP (iPopper & Etzel|[T98n) . an algorithm for analyzing 
eclipsing-binary lightcurves, and the RV model was calculated 
using our own code for solving Kepler's equation. The sizes of 
the MCMC Gaussian perturbations were set by a relatively short 
MCMC run (of ~ 10^ accepted steps). The;^^ for each trial point 
was the sum of the lightcurve x^ and the x^ of the RV data. 
Following the Metropolis-Hastings algorithm, trial points with 
lower were accepted, whereas trial points with higher were 
accepted only with probability of exip(-Ax^ /2 ), assuming the 
observational errors to be Gaussian (lFordll2005t ). 

Table |4] lists the parameters estimated by the MCMC analy- 
sis. For model parameters that showed normal posterior proba- 
bility distribution, the median value of the chain is mentioned as 
our estimation. However, for model parameters with skewed dis- 
tribution (the ones with uneven confidence limits in Table |4]) the 
most probable value is mentioned as our estimation. For most 
of the parameters, the confidence limits were estimated as the 
range of values that cover the central 68.3% of the chain. The 
lower limits of L3 and Up, however, were determined by the cho- 
sen priors, which indicates that the ground-based observations 
helped in this case for placing stronger constraints on these two 
limits. The orbital period estimated by the MCMC analysis is 
consistent with the one derived from photometric follow-up ob- 
servations (see Table [3l). 

Figure [6] presents the phase-folded RV curve of C4780 with 
our best orbital solution, and Figure |7]presents the phase-folded 
lightcurve and our best model. It can be seen that no secondary 
eclipse was detected. Given the typical noise (~ 0.25%) and the 
expected number of points inside the secondary eclipse (~ 100), 
the maximum depth of the secondary eclipse is < 0.05% (1 cr 
upper limit), which corresponds to a surface brightness ratio of 
Js < 0.06. 

Figure [8] shows the histograms of the 12 free parameters of 
the model, and Figure [9] shows some correlation plots of param- 
eter pairs that show non-zero correlation, both produced from 
the final MCMC chain. The most skewed histograms are those 
of L3, k, and r^. The most prominent (and obvious) correlation is 
between L3 and k. These facts reflect the difficulty of setting an 
upper limit for L3 directly from the photometric data, and as a 
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Table 4. The parameters of C4780 from the combined Hghtcurve and RV modeHng. 



Symbol 


Parameter name 


Value 


Units 


Orbital parameters 


P 


Orbital period 


20.6841 + 0.0006 


day 


eco^ CO 


Eccentricity x cosine longitude of periastron 


n9<i +0.002 




e sin a; 


Eccentricity X sine longitude of periastron 


-0.314 + 0.006 




K 


RV semiamplitude 


6.816 + 0.039 


kms~^ 


7 


Center-of-mass RV 


19.608 + 0.038 


kms~^ 


T 


Time of periastron* 


2454234.3 l_o.o9^°-°^ 


BJD 


f 


Mass function* 


0.000522 + 0.000010 


Mo 


Photometric parameters 


To 


Time of center of primary transit 


2454240.3144 + 0.0023 


BJD 


Js 


Surface-brightness ratio 


0.02_o.o2"'-'' 




rt 


Fractional sum of radii** (= (Ri+R2)/a) 


0.0336_o.ooo4^^-^^^^ 




k 


Ratio of radii** (=R2/Ri) 


0.095_o.ooi 




X 


Impact parameter** (= cos / • (1 - e^) • • {I + e sino;)"^) 


0.0 + 0.2 




L3 


Third-light (blending) fraction*** 


0.10_o.o3"'-'' 




Up 


Limb-darkening coefficient of primary*** 


0.57_o.oi 




i 


inclination* 


90.0 + 0.4 


degree 


Parameters estimated assuming Mi = 1.2 + 0.2 M© 




Mass ratio* (M2/M1) 


0.080 + 0.005 




a 


Semi-major axis* 


0.16 + 0.01 


AU 


M2 


Mass of secondary* 


0.096 + 0.011 


Mo 


Ri 


Radius of primary* 


1.07 + 0.07 


Ro 


R2 


Radius of secondary* 


0.104_o.oo6 


Ro 


loggi 


Surface gravity of primary* 


4.47_o.o6"'-'^ 


cgs 



* derived analytically from the fitted parameters 

** see Tamuz et al. (2006) and Mazeh et al. (2006) for the reasoning behind this definition 
*** lower limit was determined by the chosen prior (see text) 
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Fig. 6. RVs of C4780 as measured by HARPS. The solid line 
is the Keplerian model produced from the orbital parameters of 
Table IH and the dashed line is the center-of-mass velocity. The 
residuals are plotted in the lower panel. Note the diflTerent scale. 
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Fig. 7. The primary transit (left panel) and the part where the sec- 
ondary eclipse should have occurred (right panel) in the phase- 
folded Hghtcurve of C4780. Our best model is overplotted with 
a solid red line. The residuals were moved up by 0.97 for clarity. 



consequence the difficulty of setting an upper limit for the sec- 
ondary radius (even if the primary radius is known). 

The secondary mass (M2) and both radii (Ri and R2) were 
estimated using the MCMC chain and the estimated mass of the 
primary (1.2 ± 0.2 M©) in the following manner: 

- For each point of the MCMC chain the primary mass func- 
tion (/), inclination (/), and semimajor axis of the primary 
(ai) were calculated analytically. 



- A sample of 5000 Mi values, normally distributed according 
the Ml value found, was generated. 

- A subsample of 5000 points was taken from the MCMC 
chain by taking every 200th point of the chain. 

- For each Mi value of the selected sample and each of the se- 
lected MCMC points, the mass ratio (q) was calculated ana- 
lytically, using the relation 

(Mif-^ sin^ i)q^ -q^ -2q-l=0. (1) 
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Fig. 9. Selected correlation plots of parameter pairs from the final MCMC chain. For clarity, only one of every 500 points of the 
chain was plotted. 




P[day] 


vikm/s] 











0.1 0.12 0.14 0.16 0.18 0.23 0.24 0.25 0.26 20.682 20.684 20.686 19.5 19.6 19.7 19.8 




0.2 0.4 0.6 


6.7 6.8 6.9 




esinw 

4 




T -2454240 [BJD] 




_ 







0.05 0.1 0.15 0.58 0.6 0.62 0.64 -0.34 -0.32 -0.3 0.305 0.31 0.315 0.32 0.325 

Fig. 8. Histograms of the 12 free parameters of the model from 
the final MCMC chain. 



The secondary mass (M2) was then calculated for each pair 
of Ml and q values, and the semimajor axis (a) was calcu- 
lated for each set of Mi, M2, and ai values. 
Using Kepler's third law and the definitions of and k (see 
Table 131) the primary radius (Ri) can be expressed as 



^3 ^ GP\Mi + M2) 



(2) 



where G is the universal gravitational constant and P the 
orbital period. Using this expression, the primary and sec- 



ondary radii (Ri and R2) were calculated analytically for 
each Ml value of the selected sample and each of the selected 
MCMC points. Given the mass and radius of the primary, its 
surface gravity (log gi) was also calculated. 

The values and errors of all parameters calculated analyti- 
cally (/, /, a, q, M2, Ri, loggi, and R2) were estimated the same 
way as the values and errors of the free parameters of the model. 
These parameters are shown at the bottom part of Table |4] The 
highly uneven confidence limits of R2 reflect our choice of tak- 
ing the most probable value as our estimate for model parameters 
with skewed distribution. The strong correlation between k and 
L3 (Figure [9]) shows that the estimated R2 value reflects our esti- 
mation of L3 to be ~ 10%. This contamination level agrees with 
the expected one, taking the contaminants inside the CoRoT- 
photometric mask (see Section 3.1) into account. It means that 
the adopted L3 value matches the nonexistence of any other, un- 
resolved, luminous object close to C4780. This is also consistent 
with the nondetection of a secondary companion in the spectra 
with TODMOR (see Section 3.2). 



5. Discussion 



Many studies of eclipsing binaries with low-mass (0.2 - 0.8 M©) 
main- sequence stars indicated a disagreement between the- 
oretical models and observational data. In particular, it has 
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Table 5. Main properties of the stars presented in Figure [TO] 



star 


mass 


radius 


[Fe/H] 


Teff 


magnetic activity 


Porb 


V sin / 


ref. 


name 


[Mo] 


[Rq] 


[dex] 


[K] 


(Lx/Ltoi) 


[day] 


[kms-^] 




OGLE-TR123 


0.085 + 0.011 


0.133 + 0.009 


— 


— 


— 


1.8039 


— 


1 


SDSS 0857+03 


0.090 + 0.010 


0.110 + 0.004 


— 


— 


— 


0.06528 


— 


2 


J1219-39B^ 


0.091 + 0.002 


1 174+00071 


-0.21 


— 


— 


6.7600 


— 


3 


OGLE-TR122B^ 


0.092 + 0.009 


120+^-^^^ 


+0.15 


— 


— 


7.2687 


— 


4 


C4780B^ 


0.096 + 0.011 


1 n4+u.uzo 
0.149 + 0.002 


+0.20 


— 


— 


20.6841 


— 


this work 


NNSerB^ 


0.111+0.004 


— 


— 


— 


0.13008 


— 


5 


GKVirB^ 


0.116 + 0.003 


0.155 + 0.003 


— 


— 


— 


0.34433 


— 


6 


OGLE-TR106B^ 


0.116 + 0.021 


0.181+0.013 


— 


— 


— 


2.5359 


3.59 + 0.26 


7 


GJ551^ 


0.118 + 0.012 


0.141+0.007 


+0.19 


3054 


2.73E-4 + 6.5E-5 


n/a 


— 


8,9,10,11 


rl/\l- IK-ZUJ-Ui J D 


n 19/1 -1- n ni n 


U. iO/U + U.UUDU 


— u.zu 






Z.ZjU / 




1 
iz 


SDSS 0138-00 B^ 


0.132 + 0.003 


0.165 + 0.001 








0.07276 




13 


KIC1571511B^ 


141+0.005 

^•^^^-0.004 


n 17QQ+0.0014 


+0.37 






14.0225 




14 


GJ699^ 


0.146 + 0.015 


0.1869 + 0.0012 


-0.39 


3222 


2.88E-6 + 3.3E-7 


n/a 


<2.8 


8,11,15,16 


SDSS 1210+33 B^ 


0.158 + 0.006 


0.200 + 0.004 


-2.0 






0.12449 




17 


SDSS 1548+40 B^ 


0.173 + 0.027 


0.181+0.015 








0.185 




18 


RRCaeB^ 


0.1825 + 0.0131 


0.2090 + 0.0143 




3100 




0.304 




19 


2MASS 0446+19 B^ 


0.190 + 0.020 


0.210 + 0.010 




2900 




0.61879 




20 



R eferences: ( 1) Pont et al. (2006); (2) Parsons et al. (2012b); (3) Triaud et al. (20 13); (4)|P ont et al. (2005b); (5) Parsons et al. (2010); (6) 
iParsons et al. (2012c); (7) Pont et al. (2005a); (8) Bovaiian et al. (2012); (9) Demorv et al. (2009); (10) Edvardsson et al. (1993); (11) 

^JfLopez-Morales (2007); (12) Beattv et al. (2007); (13) Parsons et al. (2012a); (14) Ofir et al. (2012); (15) Lane et al. (2001); (16) 

iRoias-Avala et al. (2012) (17) Pvrzas e t all (120 12 ): (18) Pv rzas et all (2009): (19) Maxted et al. (2007); (20) Hebb et al. (2006). Notes: ^ A single 



star for which the mass was inferred from a mass-luminosity relation. ^ An sbl EB for which the mass estimation is model dependent. ^ An 
eclipsing white dwarf + M-dwarf binary for which the radius listed in the table is the volume averaged one. 




0.08- 

0.08 0.1 0.12 0.14 0.16 0.18 0.2 

Mass[M. ] 



Fig. 10. IMass-radius diagram of VLIMS (M < 0.2 Mq). The sym- 
bols represent observed stars, while the lines correspond to theo- 
retical mass-radius relations. Red rectangles are secondary stars 
of main-sequence EEs, green circles are secondary stars of white 
dwarf -h IM-dwarf EEs, and blue diamonds are single stars (all 
references are given in Table O. The black triangle stands for 
C4780E. Solid blue lines correspond to theoretical isochrones 
of solar metallicity and ages of 0.25, 1, and 5 Gyr (goin g from 
top to bottom along the 0.08 Mq line) from .Earaffe et aP (119981) . 
To illustrate the effect of metal licity on size dash ed-red lines 
show the Dartmouth isochrones (iDotter et al.ll200 8') of 1 Gyr for 
[Fe/H] = 0.2 (upper line) and [Fe/H] = -0.5 (lower line). 



cise masses and radii measurements of both components of CM 
Draconis ( Lacy 1977). However, in the domain of very low-mass 
stars (= VL MS, M < 0.2 M©), only a few studies have been pub- 
lished (e.g., iLopez-Morales 2007: Eoyajian et al 20 12), since 

only a handful of precise masses and radii have been derived in 
this domain. The derivation of the mass and radius of C4780E, 
albeit not a very accurate determination, is an opportunity to re- 
visit the VLMS mass-radius relation. 

In Figure [TOl we plotted the mass and radius of C4780E to- 
gether with other VLMS with derived masses and radii, the pa- 
rameters of which were given in Table |5] together with a few 
available VLMS models. The figure suggested that C4780E 
might be the smallest main- sequence star detected so faiQ. Its ra- 
dius is consistent and might even be below the radius predicted 
by theoretical models for an M-star with such mass, metallicity, 
and age. 

The models presented in Figure [TO] display a theoretical 
spread, probably caused mainly by diff'erences in age and metal- 
licity. Several additional parameters, not plotted in Figure [TOl 
can have further impact on the mass-radius relation. In particu- 
lar, fast rotation, magnetic activity, strong irradiation, and clouds 
were considered in t he literature (e.g., Lopez-Morales 2007; 
"Chabrier et al. 2007; Morales et al. 2010; Eouchy et alj i2QTlt 
Eurrows et al. 2011; Knigge et al. 2011). Those eff'ects can en- 
large the theoretical spread of the models even further. Given 
this relatively large range of theoretical stellar radii, no apparent 
inconsistency between the observed systems plotted in Figure[TOl 
and the available theory can be deduced at this point. 



been claimed that the stellar radii computed from models in 
this mas s domain are 5 % - 15% lower th an observed (e.g., 
Lopez-M orales & Ribas 2005t [Morales et al . 2010; Krau s et all 
20111) . This a rgument is sometim es called the M-dwarf radius 
problem (e.g.. iTriaud et al.ll2013l) . and it dates back to the pre- 



^ iBerger et"aLl (l2009h reported a radius measurements of the primary 
in the L-dwarf binary 2MASS 0746+20 using radio emission. They de- 
rived 0.078 + 0.010 Ro for this 0.085 + O OlOMc:, L-dwarf However, 
since this result is still under debate (e.g. jKonopackv et al.l l2012), we 
decided not to include 2MASS 0746+20 in the sample listed in Table|5] 
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6. Conclusions 

C4780B is a transiting very-low-mass M-dwarf, whose ra- 
dius is consistent and might even be smaller than the predic- 
tions of theoretical models. Its discovery adds an important 
piece to the puzzle that any modem stellar-evolution theory of 
VLMS, both in binaries and as single stars, would have to ac- 
count for. Further investigation of the blending fraction (L3), 
for instance with a multicolor high-precision photometric ob- 
servations, might reduce the uncertainties on its radius, thereby 
placing even stronger constraints. Since M- stars represent an 
overwhelming fraction of the gala ctic stellar population (e.g., 
IChabriejl20Q3l: iKroupa et al.ll201 Ih . understanding their proper- 
ties is important not only for advancing stellar astrophysics but 
also for other field s of science, like terrestrial exoplanet searches 
(IScalo et aLlliOOTh . 

The size of VLMS is comparable to that of giant planets, 
making such objects in binary systems one of the main false- 
alarm sources in any transiting-planet survey (e.g., Pont et al. 
l2QQ6l: ^anterne et al. 2012). As a result, many such sys- 
tems were considered as good transiting-planet candidates and 
were followed-up spectroscopically to determine their nature. 
Complete and homogeneous analysis of the syste ms that were 
detect ed by CoRoT and Kepler, similar to what ITriaud et al.l 
(l2Q13h are aiming at using the WASP data, might considerably 
enhance our understanding of the astrophysics close to the low 
end of the main sequence. 
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